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Can photo-ionisation explain the decreasing fraction of X-ray 

obscured AGN with luminosity? 
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Abstract. Chandra and XMM-Newton surveys show that the fraction of obscured AGN decreases rapidly with 
increasing luminosity. Although this is usually explained by assuming that the covering factor of the central engine 
is much smaller at luminous QSOs, the exact origin of this effect remains unknown. We perform toy simulations to 
test whether photo-ionisation of the obscuring screen in the presence of a strong radiation field can reproduce this 
effect. In particular, we create X-ray spectral simulations using a warm absorber model assuming a range of input 
column densities and ionisation parameters. We fit instead the simulated spectra with a simple cold absorption 
power-law model that is the standard practice in X-ray surveys. We find that the fraction of absorbed AGN 
should fall with luminosity as L-° i6±o.03 in h 

agreement with the observations. Furthermore, this apparent 
decrease in the obscuring material is consistent with the dependence of the FeKa narrow-line equivalent width 
on luminosity, ie. the X-ray Baldwin effect. 
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1. INTRODUCTION 

The deepest Chandra X-ray surveys find a high surface 
density of AGN (Alexander et al. 2003, Giacconi et al. 
2002). At the faintest fluxes probed, f x (2 - lOkeV) « 
10 -16 erg cm -2 s _1 , the large majority of these sources 
(~ 80%) are obscured, having a hydrogen column den- 
sity of N H > 10 22 cm- 2 (e.g. Akylas et al. 2006; Tozzi 
et al. 2006). The high fraction of obscured AGN wit- 
nessed by Chandra is well in agreement with previous 
X-ray background synthesis models (e.g. Comastri et al. 
1995, 2001). However, the situation is drastically different 
at the brighter fluxes probed by the XMM-Newton sur- 
veys. These show a marked deficit of X-ray obscured AGN 
(Piconcelli et al. 2002; Perola et al. 2004; Georgantopoulos 
et al. 2004; Caccianiga et al. 2004). The fraction of the ob- 
scured AGN drops to - 30% at fluxes of f x (2 - 10 keV) w 
10~ 14 erg cm~ 2 s _1 . In flux limited samples, a small in- 
crease in the obscured AGN fraction with decreasing ob- 
served flux is expected, as the absorbed sources present 
lower flux. However the observed increase is much steeper 
than that predicted (eg. Piconcelli et al. 2002). 

This discrepancy has been resolved with recent 
Chandra and XMM-Newton observations (e.g. Akylas et 
al. 2006; La Franca et al. 2005; Ueda et al. 2003) which 
demonstrate that the X-ray obscured AGN fraction is 
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lower at higher luminosities. This can also explain the 
deficit of type-2 narrow-line QSOs at high redshift (Steffen 
et al. 2003). Recent X-ray background synthesis models 
that take into account the dependency of the obscured 
AGN fraction on luminosity (Gilli, Comastri & Hasinger 
2007) are very successful in predicting the observed flux 
distribution. 

The important question is what is the physical mecha- 
nism that can produce the decrease in the fraction of ob- 
scured AGN. It is possible that in the presence of a strong 
radiation field the fraction covered by the torus may be 
reduced. Koningl & Kartje (1994) for example proposed a 
disc-driven hydromagnetic wind model of the torus where 
for high bolometric luminosities the radiation pressure is 
expected to flatten the torus. A similar model, the reced- 
ing torus model, has been invoked by Lawrence (1991) 
and Simpson (2005), to explain the fraction of radio AGN 
optically classified as type-2. In the presence of a strong 
radiation field the dust is heated and sublimates when 
it reaches a temperature of about 1500-2000K. At high 
luminosities the sublimation radius increases and thus ef- 
fectively the torus opening angle increases, resulting in a 
higher number of type-1 AGN. 

A similar effect is expected to take place in X-ray wave- 
lengths. The X-ray absorbing gas is situated close to the 
nucleus as inferred from X-ray variability studies of type- 
2 AGN (Risaliti, Elvis & Nicastro 2002). Thus it is rea- 
sonable to assume that a strong radiation field will also 
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ionise an appreciable fraction of the surrounding obscuring 
screen, resulting in the decrease in the "effective" column 
density. In this paper, we perform toy X-ray spectral sim- 
ulations to estimate the strength of this effect. We also 
discuss any possible connection with the X-ray Baldwin 
effect. 

2. X-ray Simulations 

2.1. Simulation Description 

We assume a constant n e r 2 relation for all AGN, where 
n e and r are the density of the obscuring screen and its 
distance from the nucleus respectively. This choice ensures 
a uniform ionisation state for the circumnuclear material 
i.e. the ionisation parameter £(=L/n e r 2 ergs cm s" 1 ) is 
the same at all distances from the nucleus. This could 
also imply that at a given distance from the central Black 
Hole the density is roughly the same in all AGN. Although 
this may not be the case, it is indirectly supported by the 
observational estimates for the density of the Broad Line 
Region (BLR) suggesting a uniform value of ~10 9 cm -3 
(Peterson 2006) across a large number of nearby AGN. 

We use the XSPEC vl2.0 package (Arnaud 1996) to 
generate 1000 simulated Chandra spectral files according 
to the following specifications: 

1. A power-law model modified by an ionised absorber 
(ABSORI model, Done et al. 1992) with a photon in- 
dex fixed to 1.9 (Nandra & Pounds 1994). 

2. A redshift of 0.8 that corresponds to the peak red- 
shift in the deep Chandra X-ray fields (e.g. Tozzi et al. 
2006) and thus the average redshift of the sources that 
produce the X-ray background. 

3. An amount of Nh randomly varying in the range of 

10 21 — 10 24 cm" 2 . We use a ratio of obscured (Nh > 

10 22 cm" 2 ) to unobscured (N H < 10 22 cm" 2 ) sources 
of 4:1 that is similar to that observed in the Chandra 
deep fields at the lowest flux levels (e.g. Akylas et al. 
2006, Tozzi et al. 2006). We further assume, a flat dis- 
tribution of N H in the 10 22 to 10 24 cm 2 region. 

4. An ionisation parameter, £=L/n e r 2 ergs cm s , 
(where L is the source luminosity in the 13.6 eV to 
13.6 keV bandpass), randomly varying between and 
1000. 

5. An arbitrary combination of an intrinsic 2-8 keV flux 
of ~ 10" 14 erg cm" 2 s" 1 and an exposure time of 50 
ks to ensure a minimum of 10 counts per spectrum 
regardless of the Nh and £ values. 

These files are then fitted in XSPEC, using a power-law 
model modified by a cold photoelectric absorption model 
( WA *PO) to estimate the "effective" neutral column den- 
sity. We consider only the 0.3-8 keV energy band in the 
spectral fits where the Chandra effective area is high. 
Sources with adequate photon statistics (> 100 counts) 
are fitted using the x 2 statistic while for sources with less 
than 100 counts we use the most appropriate C-statistic 
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Fig. 1. The solid line histogram describes the input Nh 
distribution used for the warm absorption simulations and 
the short dashed line the output one, after fitting the sim- 
ulated data with a cold absorption power-law model. 



technique (Cash 1979). We compare the parent and mea- 
sured column density distribution in Fig. [I] The solid line 
histogram describes the input Nh distribution assumed 
in the simulation (ABSORI*PO model) and the short 
dashed line the output one (WA*PO model). Note that 
due to the poissonian errors in each Nh bin the shape of 
input Nh slightly deviates from the assumed one. Clearly, 
there is a significant reduction in the number of obscured 
(Nh > 10 22 cm" 2 ) sources in comparison to the initial 
distribution. 

2.2. The simulated spectra 

An important question that arises is whether the typical 
X-ray surveys such as carried out by Chandra or XMM- 
Newton can discriminate between a photo-ionised (i.e. 
ABSORI model) and a neutral absorption model (i.e. WA 
model) given the observational limitations i.e. the limited 
signal-to-noise ratio and energy resolution. The answer to 
this question depends primarily on the spectrum quality 
in the sense that in the sources with many counts we can 
more easily reject the neutral absorption model based on 
the x 2 statistic result. The amount of the Nh and its de- 
gree of ionisation £ also play an important role given that 
their value may cause more or less prominent ionisation 
features to appear in the X-ray spectrum. 

To explore this we examine two different samples, 
each one representative of our simulation, containing ^100 
sources each. These sets present ^100 and ^500 counts re- 
spectively in their 0.3-8 keV spectrum. In Fig.[2]we present 
the distribution of the best reduced x 2 values (x 2 ) when 
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The xl distribution for sources with ~100 counts 
The xl distribution for sources with ~500 counts 



Table 1. Spectral fitting results for the eight example 
simulated spectra presenting in Fig [3] 
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Fig. 2. The distribution of the xt values for the "low 
count" (dotted line) and the "high count" (solid line) sam- 
ples. 



fitting the data using a cold absorption power-law model 
with T fixed to 1.9. 

Both distributions presented in this figure peak around 
X 2 ~ 1 while the high count distribution is shifted to 
higher xt values in respect to the low count distribution. 
Furthermore these histograms suggest that the ~85 per- 
cent of the low count and the ~75 per cent of the high 
count sources would give an acceptable statistical result 
when fitted using a simple absorbed power law model (a 
value of 1 of the xt implies that the adopted model has an 
acceptance probability of 0.4 for 20 degrees of freedom). 
Only in a few cases the cold absorption power-law model 
could be automatically rejected based on the xt value (e.g. 
xt ~ 2 implies that the cold absorption power-law model 
is rejected at ~ 3cr confidence level for 20 degrees of free- 
dom). In Fig. [3] we present some characteristic high-count 
simulated spectra. Table 1 summarises the spectral fitting 
results. 

The results presented above suggest that mildly 
ionised (£ ~ 100) column densities are well described by 
a simple power law model. This is that, given the source 
count statistic, it is not easy to discriminate between a 
neutral or an ionised model. On the other hand, when 
the column density is strongly ionised (£ ~ 1000) the re- 
sulting xt values are high and thus the power-law models 
not accepted. Note however that in the simulated spectra 
there is no sign of an Fe K a edge consistent with mildly or 
highly ionised material but only a "soft excess" like feature 
at lower energies that could be modelled using a second 
steeper power law model to improve the statistical result. 
In the case of even stronger ionised spectra (£ = 5000), 
the spectral fit with a cold absorption power-law model 
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is not accepted for those sources with the higher column 
density (N H > 5 x 10 23 cm -2 (see table 1). 

However, note that these highly ionised (and thus 
highly luminous) sources exist preferentially at higher red- 
shifts and thus the low energy absorption features would 
be shifted outside the observed band. For example the 
mean redshift of the sources with Lx>10 44 erg s _1 in the 
combined sample (X MM- Newton and Chandra) of Akylas 
et al (2006) is about 1.6. This suggests that the "soft ex- 
cess" like features seen in Fig [3] would be probably shifted 
outside the observed energy band. 

2.3. Results 

In the upper panel of Fig. 2] we plot the fraction of 
obscured sources in our simulated sample as a func- 
tion of the ionisation parameter. We fit the data with a 
model of the form of F~£ a . We find that this fraction 
drops with increasing ionisation parameter as £-0- 16 ± - 03 . 
According to the definition of the ionisation parameter £, 
the same fraction may be written in terms of luminosity 
as (L/n e r 2 )~ - 16 . A similar correlation has been suggested 
by recent studies of Chandra and XMM-Newton observa- 
tions (e.g. Akylas et al. 2006; La Franca et al. 2005; Ueda 
et al. 2003). We fit the data of Akylas et al. (2006) with 
our simulation model (L /n e r 2 )~ 0A 03 leaving only the 
normalisation (i.e. (n e r 2 ) +016 ) as a free parameter The 
fit is accepted at the 80 per cent confidence level yielding 
n e r 2 ~2x 10 41 cm~ 2 cm -1 . Individual studies of the local 
bright AGN NGC1068 (e.g. Jaffe et al. 2004) and Circinus 
(e.g. Prieto et al. 2004) show that the obscuring material 
lies outside the BLR up to a distance of a few light years. 
Our estimate for the n e r 2 predicts that at these distances 
(1 pc) the density of the material should be ~10 5 cm~ 3 . 

In Fig. 2] (lower panel) we plot the observed fraction of 
obscured AGN as a function of luminosity adapted from 
Akylas et al. (2006), together with the predictions of our 
ionisation model (L~ 016 ). Our model predicts that the 
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Fig. 3. Characteristic spectra from the high-count simulated sample. We present two different Nh values, namely 10 23 
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fraction of obscured AGN reaches the assumed maximum 
of 80 per cent at a luminosity of 4x 10 42 erg s _1 . It is worth 
noting that the exact shape of the input Nh distribution 
may affect the shape of the Obscured Fraction - £ relation. 
This is true only for the sources with medium Nh values 
ie. between 10 22 and 10 23 cm -2 . Sources with Nh < 10 22 
cm~ 2 would always remain lower and sources with Nh > 
10 23 cm -2 would always be higher than the critical value 
of 10 22 cm -2 regardless of the photo-ionisation state of 
the material. Our estimations show that if we reduce the 
number of sources with 10 22 < N H < 10 23 by 50 per cent 
(see Fig. [T]) we get a flatter Obscured Fraction - £ relation 

(F-r - 1 ). 

3. Discussion 

3.1. The warm absorber 

The results presented above demonstrate that the sim- 
plest photo-ionisation model can roughly reproduce the 
decrease in the fraction of obscured AGN at high lumi- 
nosities. If this model is true it would imply that many 
luminous QSOs would present evidence for warm absorp- 
tion. The absorption edges due to OVII and OVIII at 
0.74 and 0.87 keV respectively are the strongest signa- 
tures of ionised gas at X-ray wavelengths. The high reso- 
lution achievable with the grating spectrometers on-board 
XMM-Newton and Chandra revealed many more absorp- 
tion features as e.g. the Unresolved Transitional Array 
between 0.73 and 0.77 keV (Behar, Sako & Kahn 2001). 
These features are blue-shifted and show distinct velocity 
components (eg. Kaastra et al. 2000). This strongly sug- 
gests that the warm absorbers are associated with out- 
flows. XMM-Newton high sensitivity spectroscopy with 
RGS (Guainazzi & Bianchi 2006) places tight constrains 
on the location of the warm absorber. These appear to be 
associated with the molecular torus. Warm absorbers are 
found in about 50% of nearby Seyfert galaxies (Reynolds 
1997; Crenshaw, Kraemer & George 2003). 

The picture was unclear at high (QSO) luminosities, 
are red-shifted at low energies. In the ASCA QSO samples 
of Reeves & Turner (2000) and George et al. (2000) there 
is little evidence for ionised absorption. Because of its 
high effective area, XMM-Newton routinely produced high 
signal-to-noise spectra and changed drastically the above 
picture. Piconcelli et al. (2005) find that about 50% of the 
40 QSOs in their PG sample present warm absorbers. 

It is interesting to investigate whether there is any ev- 
idence for the presence of warm absorber features in the 
spectra of luminous AGN in deep surveys. Here we exam- 
ine the luminous AGN in the combined sample of Akylas 
et al. (2006). There are about 30 QSOs with L x > 10 44 
erg s _1 below a redshift of <~1.2 where the oxygen ab- 
sorption edges should be readily detected. Of these only 
10 have a good signal-to-noise ratio with over 200 pho- 
tons in their spectra. We fit a warm absorber model and 
find that 3 sources have an equally good fit with an ionised 
£ > 100 warm absorber, as compared to the cold absorber. 




42 43 44 45 

Log[L] (2-8 keV) ergs s" 1 

Fig. 4. Upper panel: The simulation results for the frac- 
tion of obscured AGN as a function of the ionisation 
parameter (points) together with the best fit power-law 
model, F^ £-ci6 ( so lid line). Lower panel: The frac- 
tion of the obscured AGN as a function of luminosity 
in the combined XMM-Newton and Chandra data pre- 
sented in Akylas et al. (2006) (points). The solid line rep- 
resents the correlation predicted by our simulation model 
(F~L-°- 16 ). 



Obviously the limited photon statistics do not allow for a 
conclusive test in this case. 



3.2. The X-ray Baldwin effect 

The Fe lines above 6 keV provide a census of the state of 
circumnuclear material in AGN. It is therefore reasonable 
to expect a relation between the obscuration fraction and 
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the strength of the Fe lines. Iwasawa & Taniguchi (1993) 
first find that the equivalent width (EW) of the narrow 
Fe Ka at 6.4 keV line decreases with the 2-10 keV X- 
ray luminosity in a sample of Seyfert galaxies observed by 
Ginga (the so called Iwasawa- Taniguchi or X-ray Baldwin 
effect). The correlations above closely resemble the well 
known Baldwin effect where the EW of the CIV 1550 A 
line decreases with increasing optical luminosity (Baldwin 
1977). The optical Baldwin effect has been corroborated 
by the results of Kinney et al. (1990) and more recently 
Croom et al. (2002) using 2dF data. The last authors find 
a relation with a slope of -0.128 ± 0.015. The Baldwin 
effect can be explained on the basis of either a diminishing 
covering factor or alternatively an increasing ionisation 
parameter (Mushotzky & Ferland 1984) with increasing 
luminosity. 

Page et al. (2004) again find evidence for an X-ray 
Baldwin effect using Seyfert and QSO data that span 
more than 5 orders of magnitude in luminosity. They 
find a correlation between the narrow core of the Fe K a 
EW and the 2-10 keV X-ray luminosity in the form of 
EW oc L x 017±0 08 . More specifically, the EW obtains val- 
ues around 150 eV at low luminosities (Lx(2 — lOkeV) < 
10 43 erg s^ 1 ), dropping to less than 50 eV at high lumi- 
nosities (Lx > 10 46 erg s _1 ). Jiang et al. (2006) confirmed 
the X-ray Baldwin effect using Chandra data. Recent re- 
sults by Bianchi et al. (2007) shed more light to the above 
issues. They find a highly significant anti-correlation of 
narrow Fe K a line EW with respect to the Eddington 
ratio while no dependence on the black hole mass is ap- 
parent. Moreover they find a correlation between the ratio 
of highly ionised to neutral FeK line flux with luminosity 
that supports our photo-ionisation scenario. Note however 
that according to their claims this is a rather weak result 
(Log(R) = (0.12±0.06) + (0.09±0.06)Log(Lx 44 ) to explain 
the observed X-ray Baldwin effect. 

It is interesting to note that the X-ray Baldwin effect 
is also present in the broad component of the FeK line. 
Nandra et al. (1997) find an anti-correlation of the broad 
Fe line EW with the 2-10 keV luminosity in a sample of 18 
Seyfert galaxies observed by ASCA. The broad FeK line 
most likely originates in the accretion disk. Nandra et al. 
(1997) assert that the most obvious explanation for the 
X-ray Baldwin effect observed is the photo-ionisation of 
the accretion disk. 

4. Conclusions 

We employ X-ray spectral simulations to determine the ef- 
fect of photo-ionisation on the obscuring torus in AGN. In 
particular, our goal is to explore whether photo-ionisation 
can explain the steep decrease in the fraction of obscured 
AGN with increasing luminosity as witnessed in X-ray sur- 
veys. Our conclusions can be summarised as follows: 

— Photo-ionisation can reproduce the observed decrease 
in the absorbed AGN fraction as a function of luminos- 
ity. We find that the fraction of absorbed AGN should 



fall with luminosity as £-°- 16 ±o.03 j n rough agreement 
with the observations. 

— The observed decrease in the obscuring material cov- 
ering factor is roughly consistent with the dependence 
of the Fe K a line equivalent width on luminosity, the 
X-ray Baldwin effect. 

— A simple prediction of the photo-ionised model is that 
a high fraction of luminous QSOs should present warm 
absorption features in their X-ray spectra. Future, high 
effective area mission, such as XEUS will be able to 
quantitatively assess this scenario. 
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